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Abstract — The Moderate Resolution Imaging Spectroradiome- 
ter (MODIS) instruments currently operate onboard the National 
Aeronautics and Space Administration (NASA’s) Terra and Aqua 
spacecraft, launched on December 18, 1999 and May 4, 2002, 
respectively. MODIS has 36 spectral bands, among which 20 are 
reflective solar bands (RSBs) covering a spectral range from 0.412 
to 2.13 /im. The RSBs are calibrated on orbit using a solar 
diffuser (SD) and an SD stability monitor and with additional 
measurements from lunar observations via a space view (SV) 
port. Selected pseudo-invariant desert sites are also used to track 
the RSB on-orbit gain change, particularly for short-wavelength 
bands. MODIS views the Earth surface, SV, and the onboard 
calibrators using a two-sided scan mirror. The response versus 
scan angle (RVS) of the scan mirror was characterized prior 
to launch, and its changes are tracked using observations made 
at different angles of incidence from onboard SD, lunar, and 
Earth view (EV) measurements. These observations show that 
the optical properties of the scan mirror have experienced large 
wavelength-dependent degradation in both the visible and near 
infrared spectral regions. Algorithms have been developed to track 
the on-orbit RVS change using the calibrators and the selected 
desert sites. These algorithms have been applied to both Terra 
and Aqua MODIS Level IB (LIB) to improve the EV data accu- 
racy since LIB Collection 4, refined in Collection 5, and further 
improved in the latest Collection 6 (C6). In C6, two approaches 
have been used to derive the time-dependent RVS for MODIS 
RSB. The first approach relies on data collected from sensor 
onboard calibrators and mirror side ratios from EV observations. 
The second approach uses onboard calibrators and EV response 
trending from selected desert sites. This approach is mainly used 
for the bands with much larger changes in their time-dependent 
RVS, such as the Terra MODIS bands 1-4, 8, and 9 and the Aqua 
MODIS bands 8 and 9. In this paper, the algorithms of these 
approaches are described, their performance is demonstrated, and 
their impact on LIB products is discussed. In general, the shorter 
wavelength bands have experienced a larger on-orbit RVS change, 
which, in general, are mirror side and detector dependent. The 
on-orbit RVS change due to the degradation of band 8 can be as 
large as 35% for Terra MODIS and 20% for Aqua MODIS. Vital 
to maintaining the accuracy of the MODIS LIB products is an 
accurate characterization of the on-orbit RVS change. The derived 
time-independent RVS, implemented in C6, makes an important 
improvement to the quality of the MODIS LIB products. 
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I. Introduction 

T WO Moderate Resolution Imaging Spectroradiometer 
(MODIS) instruments are currently on orbit, one aboard 
the Terra spacecraft launched on December 18, 1999 and the 
other on the Aqua spacecraft launched on May 4, 2002 [1]- 
[4]. MODIS has 36 spectral bands ranging in wavelength from 
0.412 to 14.4 /im. Continuous global observations are made 
at three nadir spatial resolutions: 0.25 km for bands 1 and 2, 
0.5 km for bands 3-7, and 1 km for the bands 8-36 [5], [6]. 
The detectors of the 36 bands are distributed on four focal 
plane assemblies (FPAs): visible (VIS), near infrared (NIR), 
short-wave infrared/midwave infrared, and long-wave infrared. 
MODIS has four onboard calibrators in addition to a space view 
(SV) port. The calibrators are a solar diffuser (SD) with an SD 
stability monitor (SDSM) [7]-[9], a blackbody (BB) [10], [11], 
and a spectroradiometric calibration assembly (SRCA) [12], 
[13]. MODIS views the Earth surface, SV, and the onboard 
calibrators via a rotating two-side scan mirror. 

The MODIS scan mirror is shown in Fig. 1. It is constructed 
from a polished nickel-plated beryllium base coated with high- 
purity silver, which is then overcoated with the Denton pro- 
prietary silicon monoxide and silicon dioxide mixture [14]. 
Using both sides of the scan mirror, MODIS collects data over 
five sectors, the SD, BB, SRCA, SV, and Earth view (EV), 
sequentially with a scan period of 1.477 s. The light from each 
of these sectors is incident on the scan mirror at a different angle 
of incidence (AOI). The AOI varies from 10.5° to 65.5° for EV 
with the SD located at 50.2°, the BB at 26.2°, the SRCA (nadir) 
at 38.2°, and the SV at 1 1.2° as shown in Fig. 2. The gain of the 
instrument at any AOI for a given band, detector, and mirror 
side (MS) depends on the reflectance of the MS at the AOI, the 
optical elements in its path, the electronic configuration, and the 
detector. Since the last three items are independent of the AOI 
of the scan mirror, the variation of the instrument gain with the 
AOI strongly depends on the reflectance of the scan mirror. A 
response versus scan angle (RVS) function is used to describe 
the dependence of the MODIS response on the AOI for each 
band, detector, and MS. 

Among the 36 bands, bands 1-19 and 26 are the reflec- 
tive solar bands (RSBs). Their wavelengths and specifications 
are listed in Table I. The RSB RVS of the scan mirror for 
both MODIS instruments was characterized prior to launch 
by measuring the response at different AOIs using a spherical 
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Fig. 2. AOIs for the MODIS scan mirror. 

integrating source (SIS) [15]— [17]. The derived prelaunch RVS 
clearly indicates a dependence on both AOI and wavelength. 
The RSBs are calibrated on orbit by the SD with its degra- 
dation being tracked by the SDSM [7]— [9] . Both instruments 
are scheduled to view the Moon approximately monthly [18], 

[19] . The lunar calibration can be used to derive the RSB on- 
orbit gain change at the AOI of the SV. MODIS has monthly 
SRCA measurements in radiometric mode, which provide some 
information for the RSB on-orbit gain change at the AOI of 
the nadir [12], [13]. The selected ocean sites and pseudo- 
invariant desert targets are also used to track the MS ratios and 
the response trending at various AOIs, respectively. The RVS 
on-orbit change can be derived using all the aforementioned 
measurements with approximation for the AOI dependence 

[20] — [24] . 

The time-dependent RSB RVS algorithms have been devel- 
oped and applied in MODIS Level IB (LIB) since Collec- 
tion 4 (C4). The C4 algorithms were developed based solely 
on the measurements of the onboard calibrators, which also 
include lunar observations, and the time-dependent RVS was 
applied to short-wavelength bands exhibiting a relative large 
change in their responses. In MODIS LIB Collection 5 (C5), 
the algorithms were refined to include the EV MS response 
ratios in the derivation of the time-dependent RVS. In the 
latest MODIS LIB version, Collection 6 (C6), the on-orbit 
RVS is characterized using a combination of measurements 
from the onboard calibrators as used in C5 and EV response 


trending at other AOIs from pseudo-invariant desert targets 

[21]. For most MODIS RSBs with relatively small changes in 
their responses, the time-dependent RVS can still be derived 
using the onboard SD and lunar observations. In the case of 
short-wavelength RSB, the SD/SDSM and lunar observations 
alone may not provide an accurate estimation of the gain, and 
furthermore, the gain at other AOIs is required to determine 
the on-orbit RVS change. For these bands, the EV response 
trending over pseudo-invariant desert targets is included to 
derive the RVS. In this paper, the formulation of the MODIS 
C6 RSB RVS algorithms is presented, and the performance of 
the on-orbit RVS algorithms is discussed. Improvements in the 
LIB products with the time-dependent RVS derived using the 
on-orbit RVS algorithm are also demonstrated. The bands with 
shorter wavelength have exhibited larger RVS change on orbit. 
For certain wavelength bands with large RVS change, the RVS 
has also exhibited detector dependence, particularly for Terra 
MODIS. In spite of being identical instruments. Terra and Aqua 
MODIS have exhibited significantly different on-orbit behavior 
of the RVS change. In general, the Terra RSBs have experienced 
more change, primarily due to the SD door anomaly (July 2, 
2003) whereby the SD door was set to a permanent open 
position. In addition, due to the smoking incident during the 
prelaunch characterization [8], the difference in the responsivity 
of the two MSs of Terra MODIS RSB has exhibited significant 
variation on orbit. 

II. RVS 

A. Prelaunch RVS 

The prelaunch RVS measurements for both Terra and Aqua 
MODIS RSBs were performed using the SIS to illuminate the 
scan mirror at 13 different AOIs with three lamp illumination 
levels, high for bands 1, 3, 4, 8, and 9, low for bands 13-16, 
and medium for all other bands. Ten scans of data with each 
containing about 100 frames were collected during each mea- 
surement [17]. The background-subtracted instrument response 
is averaged over the ten center frames and five scans with the 
same MS for each band B and detector D. Then, the averaged 
background-subtracted MODIS response to the illumination of 
the SIS for each band, detector, and MS at 13 AOIs is fitted to 
a quadratic form 

/bdm( 0) = «r M + c^ DM e + af DM 6 2 (1) 

where 9 is the AOI. The RVS for MODIS RSB is normalized 
at the AOI of its onboard SD. The derived prelaunch RVS has a 
clear band (or wavelength) and MS dependence but a negligible 
detector dependence (within 0.1%). Hence, a band-averaged 
prelaunch RVS was used for MODIS RSB, which was obtained 
using 

r BM (^) = -Jj- /bdm(0)/bdm(0sd) (2) 

D D 

where Nd is the number of detector of the band and $sd is the 
AOI oftheSD. 

Fig. 3 shows the prelaunch RVS for Terra bands 3, 8, and 17. 
A clear dependence on the AOI and MS is seen in the prelaunch 
RVS of all three bands. This indicates that the two MSs of 
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TABLE I 

MODIS RSB Key Specification and Basic Information 


Band 

CW 

BW 

IFOV 

FPA 

Ltyp 

SNR 

Primary Use 

1 

0.645 

50 

0.25 

NIR 

21.8 

128 

Land/Cloud/Aerosols Boundaries 

2 

0.858 

35 

0.25 

NIR 

24.7 

201 

3 

0.469 

20 

0.5 

VIS 

35.3 

243 


4 

0.555 

20 

0.5 

VIS 

29.0 

228 


5 

1.240 

20 

0.5 

SMIR 

5.4 

74 

Land/Cloud/Aerosols Properties 

6 

1.640 

24 

0.5 

SMIR 

7.3 

275 


7 

2.130 

50 

0.5 

SMIR 

1.0 

no 


8 

0.412 

15 

1.0 

VIS 

44.9 

880 


9 

0.443 

10 

1.0 

VIS 

41.9 

838 


10 

0.488 

10 

1.0 

VIS 

32.1 

802 


11 

0.531 

10 

1.0 

VIS 

27.9 

754 


12 

0.551 

10 

1.0 

VIS 

21.0 

750 

Ocean Color/Phytoplankton/ 
Biogeochemistry 

13 

0.667 

10 

1.0 

NIR 

9.5 

910 

14 

0.678 

10 

1.0 

NIR 

8.7 

1087 


15 

0.748 

10 

1.0 

NIR 

10.2 

586 


16 

0.869 

15 

1.0 

NIR 

6.2 

516 


17 

0.905 

30 

1.0 

NIR 

10.0 

167 


18 

0.936 

10 

1.0 

NIR 

3.6 

57 

Atmospheric Water Vapor 

19 

0.940 

50 

1.0 

NIR 

15.0 

250 


26 

1.375 

30 

1.0 

SMIR 

6.0 

150 

Cirrus Clouds Water Vapor 



Fig. 3. Prelaunch RVS for Terra bands 3, 8, and 17. 

the Terra MODIS scan mirror exhibit significantly different 
behavior prior to launch. For band 8 (0.412 /im), the prelaunch 
RVS has a strong AOI dependence with a 5% variation over the 
EV scan-angle range. In the case of bands 3 (0.469 /im) and 17 
(0.905 /im), the prelaunch RVS exhibits a change of less than 
1% over the scan-angle range. Fig. 4 shows the prelaunch RVS 
for Aqua bands 3, 8, and 17. For these bands, the prelaunch RVS 
MS differences are less than 0.1%, much smaller than those 
observed in Terra RSB RVS. Their AOI dependence is also less 
evident in comparison to that of Terra MODIS. The maximum 
variation of up to 1.5% is observed in the prelaunch RVS for 
Aqua band 8 across the entire scan-angle range. This variation 
is less than one-third of the corresponding behavior observed in 
Terra MODIS [17], 

B. RVS On-Orbit Update Algorithms 

The RVS of the scan mirror may change after prelaunch 
measurements, particularly after launch. The variations in the 
RVS from prelaunch characterization to the first on-orbit mea- 



Fig. 4. Prelaunch RVS for Aqua bands 3, 8, and 17. 

surements cannot be accurately characterized. Considering that 
the two MODIS instruments were well preserved before launch, 
the aforementioned variations of the RVS of the MODIS RSB 
are neglected in this analysis. Then, we can express the RVS as 

r-BDM(M) = 4 L(%b°dm(M) (3) 

where t) denotes the on-orbit RVS change for band B, 

detector D, and MS M from the time t o (time when the nadir 
door was first opened) at which equals 1, i.e., 

r“ DM (0, t 0 ) = 1. Since the RSB RVS is normalized at the AOI 
of the SD, the on-orbit RVS change can be expressed as 

fBDM(^f) = 3BDm(^i t)/ <?BDm($SDj t) (4) 

where <7BDM(0 5 t) anc * 5BDM(#SD,t) are the on-orbit gain 
change at the corresponding AOI and that of the SD, 
respectively. 

The MODIS RSB on-orbit gain change can be tracked by SD/ 
SDSM calibration, lunar observations, and response trending 
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over selected EV sites. The knowledge of the absolute radiance 
of the observed objects may not be necessary for the on-orbit 
gain change tracking, but their reflectance must be stable during 
the relevant period of time. The response ratios of the two MSs 
at any AOI can also be obtained from the selected EV sites, 
which do not require long-term stability of the reflectance. The 
ratios can be used to derive the MS ratios of the on-orbit gain 
change, which enables the determination of the on-orbit RVS 
change for one MS from that of the other side. As the response 
trending over EV sites is usually noisier than that obtained 
from the on-orbit calibrators, there is an additional uncertainty 
associated with it. For most MOD1S RSBs, using the onboard 
calibrators and the EV MS ratios is sufficient to track the on- 
orbit RVS change and is the preferred method. However, for 
some RSBs, particularly the short-wavelength bands, the EV 
response trending is needed to accurately track the on-orbit 
RVS change due to either large uncertainty of the SD calibration 
or the necessity to track the on-orbit gain change at multiple 
AOIs or both. In the following sections, two approaches are 
presented. One is based on the measurements of the onboard 
calibrators and the EV MS ratios, and the other is based on the 
lunar and the response trending over the pseudo-invariant sites. 
The former is referred to as Approach 1, and the latter is referred 
to as Approach II. 

1) On-Orbit RVS Change Derivation Using Onboard Cali- 
brators and EV MS Ratios (Approach I): As discussed earlier, 
the MODIS RSB gain at the AOI of the SD is tracked using 
the SD calibration [7]— [9] . By normalizing the gain obtained 
at the time that the nadir door was first opened, the on-orbit 
gain change at the AOI of the SD, (?bdm ($sd > i), can be 
derived from the SD calibration. The lunar calibration that 
occurs approximately monthly can provide the RSB on-orbit 
gain change, i/bdm($sVj £), at the AOI of the SV. Then, the 
RVS at the AOI of the SV can be calculated using 

r BDM(^SV,i) = ffBDM^SV^VffBDM^SDjf)- (5) 


Normalizing the RVS to the AOI of the SD, the on-orbit RVS 
change at any AOI, r\'» m ((). : t), can be expressed using a linear 
approximation for the AOI dependence as follows: 

r B°DM(M) = 1 + 77 [^BDM^SV, t) — 1] ■ (6) 

VSV ~ y SD 


The instrument EV response MS ratios can be derived from 
ocean as well as desert sites as long as the targeted areas are 
spatially uniform during each observation. Compared to ocean 
sites, a desert target yields a higher signal-to-noise ratio (SNR) 
and is deemed more suitable. However, most high-gain ocean 
bands (bands 8-16) saturate over desert surfaces. Consequently, 
scenes from oceans are chosen to calculate the MS ratio [20]. 
Using an ocean target, the instrument response MS ratios can 
be obtained at multiple AOIs and can be expressed as 


SBd(®i, t) 


dn{B,D,M2,6j,t) . 

dn(B,D,Ml,Oi,ty * 1 


(7) 


where Ml and M2 denote the two sides and N\ is the number 
of Oi. 


With the MS ratios calculated, the RVS on-orbit variation of 
the two MSs can be expressed as 


r °° 

' BDM2 


(6M) 


la ^2BDMl(0SD,f)rg M i(6>i,f) M oo {D ^ ,<n 

— SBD(Pi,r) — — j — -tBDMl(.Pii U- (°) 

fl , BDM2(C'SD, t)r BM2 (Vi, U 


Thus, if the RVS on-orbit change for Ml is known, the 
variation for M2 can be calculated at multiple AOIs using a 
suitable polynomial approximation. The linear approximation 
for the AOI dependence of i'bdmi(^) ma y not adequate 
as ?'bd M 2 (^;^) a l so exhibits noticeable nonlinear dependence 
for most RSBs. The nonlinear behavior of the two MSs may 
not be identical, which may induce MS-based striping in the 
calibrated EV images. With (8) applied to calculate the RVS 
on-orbit variation for M2, the two MSs are calibrated to about 
the same uncertainty, and the striping effect due to the MS 
difference is reduced. For both Terra and Aqua MODIS RSBs, 
the MS ratios are calculated at ten AOIs and fitted to a quadratic 
form with a least-mean-square-fit approach, which is then used 
to derive the on-orbit RVS change for MS2 [21], [22]. 


2 ) On-Orbit RVS Change Derivation Using Desert Trending 
and Lunar Calibration (Approach II): Since both MODIS in- 
struments have been operating beyond their designed lifetime of 
six years, there have been concerns over the accuracy of the SD- 
based calibration and the validity of the linear approximation 
for on-orbit RVS characterization, although the SD calibration 
and the RVS linear approximation have worked reasonably well 
for the first few years. To tackle these concerns, an alternative 
source in the form of pseudo-invariant desert sites [21] has 
been used to trend the MODIS RSB response. These sites are 
considered to be temporally stable and are endorsed by the 
Committee on Earth Observation Satellites (CEOS) for sensor 
calibration/intercomparison. In principle, one site can be used 
to track the MODIS RSB on-orbit change at 16 different AOIs 
since the MODIS orbit has 16-day repeating cycles. However, 
not every measurement is useful since the reflected sunlight can 
be impacted due to the presence of clouds above the site during 
observation. Therefore, multiple sites are needed to track the 
on-orbit gain changes at various AOIs. Three widely used sites, 
Libya 1 (24.42 N, 13.35 E), Libya 2 (25.05 N, 20.48 E), and 
Libya 4 (28.55 N, 23.39 E), are selected to track the RSB on- 
orbit gain changes for both Terra and Aqua MODIS. With these 
three sites, the on-orbit gain change can be effectively tracked 
at 12 AOIs for both instruments. Since the two instruments have 
different orbits, the 12 AOIs are not identical, but they are close. 

As mentioned earlier, lunar calibration can track the RSB 
on-orbit gain changes at the AOI of the SV. Because the 
lunar surface is very stable, the lunar calibration should be 
reliable and accurate. An exception is the short-wavelength 
infrared (SWIR) bands which have strong crosstalk among each 
other. For a general discussion, assume that the on-orbit gain 
changes are derived from the selected desert sites and lunar 
calibration at N% different AOIs and use the terms .r/ r j \r ( D i , t). 
ffBM(02,i), • • • , <?BM(0JV 2 )i) t0 denote the on-orbit gain 
changes. Since the desert trending produces noisier results than 
that obtained from the onboard calibrators, only the detector- 
averaged response is derived from the EV measurement. For 
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both Terra and Aqua MODIS, N 2 is 12 or 13, depending 
on whether the lunar trending is included. The gains are all 
normalized to the time when the nadir door was first opened. 
At any given time t, there are N 2 gains at different AOIs. 
They can be fitted to a polynomial 

p 

h BM (e i t) = Y / b l (B,M,t)9 l (9) 

i=0 


dependence. The exact reason why the lunar trending does not 
match the desert trending at the SV AOI for Terra bands 1-2 is 
still unknown. 

With the prelaunch RVS rg M (0) and the RVS on-orbit vari- 
ation the time-dependent RVS Tbdm^j t) can be 

easily calculated for each band, MS, detector, and AOI. The 
EV frame number F, ranging from 0 to 1353, and the AOI 9 
are related by 


where P is the order of the polynomial. If the lunar measure- 
ment is used in the fitting, the polynomial is constrained to pass 
through the lunar response trending. 

Using the fitted polynomial, a detector-averaged RVS on- 
orbit variation can be derived using 


r 


OO 

BM 


(M 


^bm(^ t) 
^bm($sd, t) 


( 10 ) 


9 = 


65.5° - 10.5° 
1353 


U + 10.5°. 


( 11 ) 


For MODIS LIB, the RVS is refitted as a polynomial function 
of the EV frame number 


Q 


i= 0 


( 12 ) 


where the normalization factor /ibm(0sd,/) provides the on- 
orbit gain change at the AOI of the SD. The RVS on-orbit 
detector-dependent differences can be derived from the SD 
and lunar calibration with a linear approximation for the AOI 
dependence. 

C. RVSLUTs 

MODIS bands 5-7 and 26 are SWIR bands. No evidence 
has demonstrated that the RVS for these bands has noticeable 
change for both instruments since launch. Thus, prelaunch RVS 
has been applied to them in all MODIS LIB versions. For Terra 
bands 10-19 and Aqua bands 1-4 and 10-19, SD and lunar 
calibration sufficiently track the orbit gain change at the AOI 
of the SD and that of the SV, and the linear approximation for 
the AOI dependence of the on-orbit RVS change works well. 
For these bands, the RVS on-orbit changes are derived using 
Approach I with the MS 1 on-orbit RVS change calculated 
using (6) and that for MS 2 calculated using (8). In principle, (6) 
can be applied to either MSs. However, there was a “smoking” 
event for Terra MODIS during the prelaunch measurements 
[8], which provides some justification for the aforementioned 
choice for Terra MODIS. To be consistent, the same choice is 
applied to Aqua MODIS [22], [23]. 

For Terra MODIS bands 1M-, 8, and 9 and Aqua MODIS 
bands 8 and 9, the response trending over the pseudo-invariant 
desert sites reveals that the linearity assumption about the AOI 
dependence of the RVS on-orbit variation does not hold [2]. 
Comparison between the desert response trending around the 
AOI of the SD and the SD response trending also reveals that 
the SD measurements cannot reliably capture the on-orbit gain 
change for Terra bands 1-4, 8, and 9, particularly for the past 
few years [22]. For the aforementioned bands, Approach II is 
applied to derive the on-orbit RVS change. For these bands, the 
lunar response trending matches the desert response trending at 
the SV AOI quite well, and a quartic form is applied to describe 
the AOI dependence except for Terra bands 1 and 2. For these 
two bands, the lunar response trending does not match the 
desert response trending at the SV AOI, and only the response 
trending over desert sites is used to derive the on-orbit change 
with the same quadratic form being applied to describe the AOI 


where Q is the order of the polynomial and is band dependent. 
For Terra MODIS bands 3, 4, 8, and 9 and Aqua bands 8 and 9, 
Q is 4, while for all other bands, it is 2. The fitted coefficients of 
all the RSBs and the two MSs form a MODIS RSB RVS lookup 
table (LUT) with a time stamp t. The LUT is generated at a set 
of selected time stamps and used in MODIS LIB as input. 

Table II is a summary of MODIS RSB RVS update algo- 
rithms. The first column in the table is the band number. The 
second and sixth columns provide the methods used to derive 
the RVS: Prl denotes prelaunch RVS, while App. I and App. II 
indicate an RVS update using Approach I and Approach II, 
respectively. The third and seventh columns indicate the order 
of the polynomial used to describe the AOI dependence of the 
RVS on-orbit change at any given time. The fourth and eighth 
columns are the polynomial order used to describe the AOI 
dependence of the time-dependent RVS at any specified time. 
The fifth and tenth columns provide the information for the data 
used to update the RVS on orbit: DT for pseudoinvariant desert 
data, Moon for lunar data, SD for SD data, and OM for MS 
ratios obtained from ocean sites. 

III. On-Orbit Gain Changes 

A. On-Orbit Gain Change From Onboard Calibrators 
and EV MS Ratios 

1) SD/SDSM Calibration: During the first year on orbit, 
both Terra and Aqua MODIS performed the SD calibration on 
a weekly basis. The frequency of SD calibration was reduced 
to biweekly after the first year. Since luly 2, 2003, the SD 
door of Terra MODIS has been kept in the “open” position, 
and the SD screen is in place on a permanent basis. Since the 
SD door is always open, an SD calibration is performed every 
orbit. Prior to July 2, 2003, the SDSM was operated at the 
same frequency as the SD, and thereafter, it is operated once 
every three weeks. In the case of Aqua MODIS, the SD/SDSM 
scheduling frequency was further reduced to once every three 
weeks to prolong the lifetime of the SD/SDSM calibration 
capability. With the onboard SD/SDSM calibration, the gain 
of the MODIS electronic and optical system at the AOI of the 
SD can be derived. For each calibration event at time t, the 


3164 


IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 52, NO. 6, JUNE 2014 


TABLE II 

Information for MODIS RSB Time-Dependent RVS. DT: Desert Targets. Prl: Prelaunch. OM: Ocean Site MS Ratios 


Band 

Terra 

Aqua 

App. 

P 

Q 

Data 

App. 

P 

Q 

Data 

1 

App. II 

2 

2 

DT 

App. 1 

2 

2 

SD, Moon, OM 

2 

App. II 

2 

2 

DT 

App. 1 

2 

2 

SD, Moon, OM 

3 

App. II 

4 

4 

DT, Moon 

App. 1 

4 

4 

SD, Moon, OM 

4 

App. II 

4 

4 

DT, Moon 

App. 1 

4 

4 

SD, Moon, OM 

5 

Prl 

0 

2 


Prl 

0 

2 


6 

Prl 

0 

2 


Prl 

0 

2 


7 

Prl 

0 

2 


Prl 

0 

2 


8 

App. II 

4 

4 

DT, Moon 

App. II 

4 

4 

DT, Moon 

9 

App. II 

4 

4 

DT, Moon 

App. II 

4 

4 

DT, Moon 

10 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

11 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

12 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

13 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

14 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

15 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

16 

App. 1 

1 

2 

SD, Moon, OM 

App. 1 

1 

2 

SD, Moon, OM 

17 

App. 1 

1 

2 

SD, Moon 

App. 1 

1 

2 

SD, Moon 

18 

App. 1 

1 

2 

SD, Moon 

App. 1 

1 

2 

SD, Moon 

19 

App. 1 

1 

2 

SD, Moon 

App. 1 

1 

2 

SD, Moon 

26 

Prl 

0 

2 


Prl 

0 

2 
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Fig. 5. Gain trending for Terra MODIS bands 10, 13, and 17, derived from 
the SD calibration. 

calibration coefficient denoted as mi in this paper and in the 
literature [7] — [9] is derived for each band, detector, subframe, 
and MS. The gain at the AOI of the SD is inversely proportional 
to the calibration coefficient, and then the on-orbit gain change 
can be calculated by inversing the calibration coefficient. 

Fig. 5 shows the SD gain trending for the selected RSBs, 
which have been averaged over detectors and normalized to the 
time when the nadir door first opened. The filled symbols de- 
note the data for MS 1 while unfilled symbols are for MS 2. The 
discontinuities around day 305 (October 31, 2000) and day 549 
(July 2, 2001) are due to the changes of the MODIS electronic 
configuration from side A to side B and back to side A from side B, 
respectively. The increased rate of degradation beyond day 
1279 (July 2, 2003) is due to an increased exposure to sunlight 
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Fig. 6. Gain trending for Aqua MODIS bands 10, 13, and 17, derived from 
the SD calibration. 

after the SD door for Terra MODIS was permanently left in the 
open position. Even though the on-orbit change in the SD is 
characterized using the SDSM, the increased degradation rate 
of the optics, including the scan mirror, results in a larger gain 
change after day 1279. The gain trending for band 10 (0.488 
(im) shows more change compared to the longer wavelength 
band 17 (0.905 /in i), while band 13 (0.667 /im) shows a gain 
increase over the mission. Among the bands shown in Fig. 5, 
band 10 (0.488 /xm) has the shortest wavelength and shows the 
most change, about 30% for MS 1 and 36% for MS 2. 

The Aqua MODIS RSB gain derived from the SD is shown 
in Fig. 6 for bands 10 (0.488 /xm), 13 (0.667 /xm), and 17 
(0.905 /xm). Similar to Fig. 5, the gain coefficients have been 
averaged over detectors, fitted to piecewise smooth functions. 
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Fig. 7. Gain trending for Terra MODIS bands 10, 13, and 17, derived from 
the lunar calibration. 

and normalized to the time when the nadir door was first 
opened. Similar to the Terra band 13 (0.667 /im) trending 
shown in Fig. 5, the gain of Aqua band 13 (0.667 /im) derived 
from SD calibration also increases over time. The largest gain 
change among the three bands is seen for band 10 (0.488 /im), 
consistent with the trend observed in Terra MODIS. The gain 
for Aqua band 10 (0.488 /im) shows a decrease of about 11% 
which is of a smaller magnitude in comparison to Terra band 
10 (0.488 /im). The MS differences in the gain trending are 
also much smaller than those observed in Terra RSB. 

2) Lunar Calibration: Both Terra and Aqua MODIS have 
been scheduled to view the Moon approximately monthly, 
about nine months every year, through their SV [18], [19]. The 
lunar irradiance depends on lunar view geometry. To reduce 
the uncertainty of the view geometry effect correction, the 
lunar phase angle, defined as negative for a waxing moon 
and positive for a waning moon, is usually kept in the range 
[-56°, — 55°] for Aqua MODIS and in the range [55°, 56°] for 
Terra MODIS, respectively. Lunar calibration can be used to 
track the on-orbit gain change at the AOI of the SV as long as 
the data are not saturated during a lunar observation for each 
band, detector, and MS. 

The detector-averaged on-orbit gain changes derived from 
lunar calibration for Terra MODIS bands 10, 13, and 17 are 
shown in Fig. 7. Similar to Fig. 5, the discontinuities around 
day 305 (October 31, 2000) and day 549 (July 2, 2001) are due 
to the previously mentioned changes of the MODIS electronic 
configuration. Gain degradation rate changes around day 1279 
are due to the SD door remaining open starting on July 2, 
2003. A similar wavelength-dependent degradation pattern is 
observed from the lunar gains. Band 10 (0.488 /im) shows a 
degradation of about 41.5% for MS 1 and 42.5% for MS 2 at 
the AOI of the SV, which is about 11.5% and 6.5% more that 
those at the AOI of the SD as demonstrated in Fig. 5 for MSs 
1 and 2, respectively. Differences between the SD and lunar 
trending are also seen for other bands, particularly for short- 
wavelength bands, and justify the need for an on-orbit RVS 
characterization. 

The Aqua MODIS RSB detector-averaged gain derived from 
lunar calibration is shown in Fig. 8 for bands 10, 13, and 17. 



Fig. 8. Gain trending for Aqua MODIS bands 10, 13, and 17, derived from 
the lunar calibration. 



Fig. 9. Terra band 10 EV MS ratio trending observed at ocean sites in Pacific 
Ocean around the equator. 

Similar to the Terra results, the gain exhibits wavelength de- 
pendence since band 13 (0.667 /im) is saturated while viewing 
the Moon. Hence, a ratio approach using an unsaturated band as 
reference is applied to derive the lunar response trending [20], 
which may increase the uncertainty for the derived gains seen 
in Fig. 8. For band 10 (0.488 /im), the gain at the AOI of the SV 
decreases by about 12% and 13% for MSs 1 and 2, respectively. 
This is about 1% and 2% more than the corresponding results 
using the SD calibration method. Compared to Terra band 10 
(0.488 /im), the gain change as well as the difference between 
the two sources is much smaller. 

3) MS Ratios Over Ocean Sites: Most MODIS high-gain 
ocean color bands (8-16) saturate while observing land targets. 
Hence, ocean sites are selected to track the MODIS RSB MS 
ratios. In addition, a uniform ocean scene can provide mea- 
surements at multiple AOIs within a single granule. To reduce 
the seasonal oscillation induced by the polarization effect, the 
equator area in the Pacific Ocean with a longitude of —136° 
was used to obtain the response ratios used in calculation of 
the MODIS C6 LUT. From an EV granule which covers the 
equator area in the Pacific Ocean, the MS ratios are calculated 
at selected AOIs. 

Fig. 9 shows the trending of the EV MS ratios for Terra 
MODIS band 10 (0.48 /im) at six specific AOI. It is clearly 
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Fig. 10. Aqua band 10 EV MS ratio trending observed at ocean sites in Pacific Rg n EV response for Terra M0D IS band 8 MS 1 over the pseudoinvariant 
Ocean around the equator. desert sites without BRDF correction. 


seen that the MS ratios change with time and are strongly AOI 
dependent. A seasonal oscillation in recent years is also seen 
in the MS ratios at large AOI. This yearly oscillation in recent 
years (since 2007) is primarily due to the fact that MODIS is a 
polarization-sensitive instrument and the sensitivity is MS and 
time dependent, particularly at large AOI [27], Fig. 10 shows 
the corresponding Aqua band 10 MS ratio trending at the same 
six AOIs. In comparison with Terra band 10, Aqua band 10 
exhibits a much more stable MS response ratio. 

B. On-Orbit Gain Change From Pseudo-invariant Desert Sites 

1 ) Response Trending Over Desert Sites: The EV responses 
at different AOIs are derived from the noncalibrated Level 1A 
(L1A) covering the selected pseudoinvariant desert targets. The 
EV responses are averaged over a 20 x 20 km area centered 
at each site. The footprint of the selected 20 x 20 km region 
at 1-km nadir resolution covers one scan for each MS with ten 
along-track detectors and 20 frames in each scan. The data from 
high-resolution bands (bands 1-7) are aggregated to the 1-km 
resolution. The LI A granules with a 16-day repeat cycle are 
used to group the collected EV responses from a given site, so 
the responses in each group have a roughly identical AOI over 
the entire mission. 

The EV data are corrected for the background signal, the 
instrument temperature effects, the cosine effect of the incident 
sunlight, and the Earth-Sun distance variation over the orbits. 
The saturated pixels contaminated due to clouds are excluded. 
A uniformity filter using the combined standard errors of the 
reflectance of one VIS band (band 1, 0.645 pm) and one NIR 
band (band 2, 0.858 pm) is also applied to remove the cloud- 
contaminated data. All remaining EV pixels after the uniformity 
filtering are then averaged over detectors and 20 frames in each 
scan to derive the MS-dependent EV response. 

Fig. 1 1 shows the derived EV response trending at three dif- 
ferent scan angles (AOI) for Terra MODIS band 8 (0.412 pm) 
MS 1 . The yearly oscillation in the EV response at frame 111 
is largely the effect of the test site bidirectional reflectance 
distribution function (BRDF), which is dependent on the sea- 
sonal solar viewing geometry changes. The large slope and 
discontinuity in the early part of the mission is due to the 


change in the electronic configuration for Terra MODIS. The 
discontinuity will be eliminated after the application of the 
calibration coefficients. A smaller yearly oscillation is observed 
in the trending at the AOI of 36.6°, also affected by the site 
BRDF. There is seasonal oscillation in the trending for the 
AOI of 52.8° as well, with the amplitude of oscillation greatly 
increasing after 2008. The amplitude increase is due to the 
polarization property change of the band. In order to have a 
better assessment of the actual instrument gain change, the 
impacts of these yearly oscillations must be mitigated before 
further processing. 

A semiempirical BRDF model developed by Roujean et al. 
[25] is chosen to perform the BRDF normalization. The kernel- 
driven BRDF model is expressed as follows: 

p(A, 0, ip, ijj) = fc 0 (A) + k 1 (X)f 1 (e,(p,ip) + fc 2 (A)/ 2 (6»,^,V’) 

(13) 

where A is the band wavelength, 9, ip, and ip are the solar 
zenith angle, sensor zenith angle, and relative azimuth angle, 
respectively, and p is the modeled top-of-atmospheric (TOA) 
reflectance. The kernel fi accounts for the volume scattering, 
and the kernel / 2 represents the surface scattering and geomet- 
ric shadow casting. The site-specific BRDF coefficients kg, ki, 
and lip quantify the contributions to the TOA reflectance by 
combining the isotropic component with the two kernels de- 
scribing anisotropic components. They are related to the phys- 
ical structure and optical properties of the EV target site. The 
site- and AOI-specific BRDF coefficients can be derived from 
the linear regression to the observed MODIS TOA reflectance. 
The onboard calibration is known to be able to provide stable 
trends for the LIB product (i.e., radiance or reflectance) during 
the initial years after launch. Therefore, the observed TOA 
reflectance derived using the onboard calibrations in the first 
three years of the mission for both Terra and Aqua MODIS 
is used to determine the BRDF coefficients at each selected 
desert site and AOI. Since uncertainties are significantly larger 
to derive the coefficients that can cover the entire range of 
the MODIS scan angle due to the impact of the underlying 
atmosphere, the BRDF coefficients are actually derived for 
each individually fixed AOI using the corresponding EV data. 
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Fig. 12. EV response for Terra MODIS band 8 MS 1 over the pseudoinvariant Fig. 14. EV and lunar response for Terra MODIS band 8 MS 1 . 
desert sites after BRDF correction. 



Epoch 2000 


Fig. 13. EV response for Band 8 MS 1 at multiple sites for each AOI. 


Fig. 12 shows the Terra MODIS band 8 MS 1 EV response 
trending after the BRDF correction. Significant improvement 
in reducing the yearly oscillation is seen when compared with 
the EV responses in Fig. 1 1 . The seasonal variation may still 
be seen at large frame numbers after 2008 for Terra short-wave 
bands, as shown in Fig. 12 for band 8. The effect is primarily 
due to the polarization property change of the MODIS scan 
mirror [26], [27]. 

2) Comparison Among Different Sites and With the Lunar 
Observations: The radiometric stability of the selected sites 
needs to be investigated before the instrument responses mea- 
sured at the sites are applied to track the MODIS RSB on- 
orbit gain change. Fig. 13 shows the Terra band 8 (0.412 pm) 
responses observed at sites Libya 1 (24.42°, 13.35°), Libya 2 
(25.05°, 20.48°), Libya 4 (28.55°, 23.39°), and Mali 1 (19.12°, 
-4.85°) for AOIs of 12.0°, 12.8°, 30.8°, and 31.5°. Due to 
the wide MODIS swath and the proximity of the location of 
these desert sites, it is difficult to obtain the response from 
two different sites at the exact same AOI. In this case, the 
responses from AOIs of 12.0° and 30.8° are compared with the 
responses derived from AOIs of 12.8° and 31.5°, respectively. 
The responses of the first three years have been fitted to a 
smooth function of time and then normalized to the time when 


the nadir door was first opened. The first sets of data at AOIs 
of 12.0° and 12.8° provide similar response trending results for 
Terra band 8 (0.412 pm) at the AOI of 12.4°. The long-term 
difference between the two trending results is less than 1%. The 
last sets of data at AOIs of 30.8° and 31.5° provide the response 
trending results for Terra band 8(0.412 pm) at the AOI of 3 1 . 1 ° . 
Once again, the difference between the two trending results is 
less than 1%. Considering the uncertainty of the EV data and 
the small difference of the AOI, the two sets of trending results 
agree reasonably well. A similar conclusion can be obtained for 
Terra bands 1-4, 8, and 9 and Aqua bands 8 and 9, which also 
use the instrument responses over the desert sites to track the 
gain on-orbit changes at the selected AOI. 

It is well known that the Moon surface reflectance is very 
stable in the VIS and NIR band range. Thus, another direct 
validation of the radiometric stability of the desert sites is to 
compare the response trending observed at the desert sites with 
that obtained from the lunar observations. Fig. 14 shows the 
response change from the lunar view and desert site observation 
at the AOI of 12.0° which is very close to 1 1.2°, the AOI of the 
SV or the lunar view. The symbols are for measured data. The 
data are normalized at the time that the nadir first opened by 
applying the fitted smooth function for each of the observations. 
From Fig. 14, it can be seen that the response changes observed 
from the lunar view and desert site agree very well (within 1%). 
It is a direct validation for using the response change observed 
from the desert sites for tracking the on-orbit gain change for 
these bands. A similar trend is observed for most RSBs with 
the exception of the high-gain ocean bands and Terra bands 
1 and 2. For the high-gain ocean bands, they are saturated 
when they view the desert sites, and hence their on-orbit gain 
variation change cannot be tracked. In the case of Terra bands 
1 and 2, the EV and lunar response changes show a measurable 
disagreement in their responses. However, a contrary result 
is observed in the case of Aqua MODIS where a reasonable 
agreement (within 0.5%) is observed between the EV and lunar 
measurements. The cause for the difference between the lunar 
and EV responses of Terra bands 1 and 2 remains unknown. 
Due to the disagreement between the lunar and EV response 
trending, the lunar data are not used in the derivation of the on- 
orbit RVS change for these two bands. 
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Fig. 15. Corrected EV response for Terra MODIS band 8 MS 1 along with 
the two-piece fitting model. 


3) Measured Gain Fitted to Smooth Functions: The rate 
of the RSB gain change varies as a function of AOI, time, 
wavelength, and MS. Several smoothly connected functions 
have to be used to describe the SD and lunar on-orbit gain 
variation for each RSB. The EV response trending post-BRDF 
normalization has a large uncertainty, and the direct fitting of 
these data to piecewise smooth functions may prove inaccurate 
in describing the EV response. In order to avoid the larger 
uncertainties in the LUT, an alternative way to describe the EV 
response trending is described using via the function: 


Sbm($> t) 


(50.25° -0) 
50.25° — 11.2° 




( 0 - 11 . 2 °) 
50.25° — 11.2° 


0bm M 


(14) 


where 3BM° n (^) an d 0 bm W are the fitted smooth functions for 
lunar and SD gains normalized at the time to of the first opening 
of the nadir door. The values 50.25° and 11.25° correspond to 
the SD AOI and lunar AOI, respectively. g EM (Q, t) denotes the 
calibration coefficients derived from the SD calibration over 
Approach I RVS (mi/RVS) on-orbit change. Then, a corrected 
EV response can be defined 


^ n Ev(-®i 0i,t ) — dn EV (B, M, Ot, t ) / 0 (15) 


for band B and MS M at AOI Qi , which has the same long-term 
drift as the EV reflectance if the SD calibration coefficients 
and Approach I RVS are applied in MODIS LIB. Fig. 15 
shows the normalized dn' EW trending for Terra MODIS band 8 
(0.412 /im) at three AOIs, which are close to SV, nadir, and SD 
AOI. Compared to the Terra band 8 EV response trending in 
Fig. 12, the dn E y trending has a much simpler dependence of 
time and can be well described by a two-piece linear function, 
which is also used to fit the normalized dn EV trending for 
Terra bands 3 and 9. In the case of Terra bands 1, 2, and 
4, the corrected EV response is described using a one-piece 
linear function, whereas a one-piece quartic model is applied 
for Aqua MODIS. If represents the fitted smooth 

function for band B and MS M at AOI 0i, then the gain on-orbit 
change at AOI Qi derived from EV response can be expressed as 


(16) 


Each MODIS can view a selected desert site at multiple 
AOIs due to the 16-day orbit repeat cycle. Because the 
AOIs of different sites are usually slightly different due to 
their geolocation differences, a set of approximately evenly 
distributed AOIs for each instrument can be selected to track 
the on-orbit gain changes. 

IV. RVS and Its On-Orbit Performance 
A. On-Orbit Time-Dependent RVS 

For Terra bands 10-19 and Aqua bands 1-4 and 10-19, the 
RVS on-orbit change is derived using the SD/SDSM and lunar 
calibrations with a linear assumption for the AOI dependence 
(Approach 1). For Terra bands 1—4, 8, and 9 and Aqua bands 
8 and 9, the RVS on-orbit change is derived from the EV 
response trending over pseudo-invariant desert sites at selected 
AOI and lunar response trending with higher order polynomial 
approximation for the AOI dependence (Approach II). The re- 
sponse trending is obtained either from the onboard calibrators 
or desert sites, fitted to a piecewise smooth connected function, 
and normalized to the time that the nadir door first opened for 
each AOI, band, and MS. At any given time, the normalized 
response is used to describe the behavior of AOI dependence 
using a polynomial approximation. Then, the derived polyno- 
mial is normalized at the AOI of the SD. This normalized 
polynomial is the on-orbit RVS change at the time since the 
nadir door opened. 

The inverse of the normalization factor provides the on-orbit 
change of the calibration coefficient m i. For the Approach I 
bands, the on-orbit gain change at the AOI of the SD is also 
provided by the SD/SDSM calibration, and thus, the inverse 
of the normalization factor matches the on-orbit change of the 
calibration coefficient derived from the SD/SDSM calibration 
[7]— [9] . For bands calibrated with Approach II, the on-orbit 
gain changes are mainly tracked by the response trending over 
the pseudoinvariant desert sites, and the on-orbit change of 
the calibration coefficient provided by the normalization factor 
does not necessarily match that derived from the SD/SDSM 
calibration. The difference between them shows the inadequacy 
and uncertainty of the SD/SDSM calibration. 

Fig. 16 shows the RVS on-orbit change for Terra Band 8 
(0.412 pm) as a function of AOI at particular time stamps. 
It is clearly seen that a quartic form is needed to accurately 
describe the AOI dependence of the on-orbit RVS change. The 
solid symbols denote the on-orbit lunar measurements which 
also exhibit a good agreement with the EV response around the 
same AOI. At any given time, the factor used to normalize each 
fitted polynomial is the degradation of the gain at the AOI of the 
SD. Hence, the gain at the AOI of the SD is also determined in 
addition to the RVS using the lunar and EV response trending. 
The relative RVS curve shows the most temporal change at 
the AOI of the SV. This is also true for other bands using a 
quartic form to describe the AOI dependence. For Terra bands 
1-2, the lunar view response trending data do not match the EV 
response trending results and hence are not included in the RVS 
derivation for these bands as mentioned in Section II. For these 
two bands, a quadratic form is sufficient to fit the EV response 
trending results at the selected AOI for any given time. 


ffBM (Qi,t) — <7bm(^; *)■ 
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Fig. 16. Relative RVS for Terra MODIS band 8 MS 1. 



Fig. 17. RVS for Terra MODIS bands 10, 13, and 17 at the AOI of the SV. 

The detector-averaged RVS for MS 1 of Terra MODIS bands 
10, 13, and 17 at the AOI of the S V, generated using Approach I, 
is shown in Fig. 17. It can be clearly seen that the on- 
orbit change in the RVS for the shorter wavelength band 10 
(0.488 fin i) is about 17%, whereas the RVS changes for bands 
13 (0.667 /iin) and 17 (0.905 fin i) are about 3% and less than 
1%, respectively. The largest MS differences of the RVS are 
about 8% for band 10 (0.488 fim) and less than 1% for the 
other two bands. The Aqua MODIS RVS trends for bands 10, 
13, and 17 shown in Fig. 18 also show the shorter wavelength 
band 10 exhibiting a larger change compared to bands 13 and 17 
while the MS differences for all three bands are much smaller 
compared to the corresponding Terra bands. For Aqua bands 13 
and 17, the RVS at the AOI of the SV has generally increased 
with time, in contrast with Aqua band 10 (0.488 /iin) and Terra 
bands 10, 13, and 17. As mentioned in Section II, the RVS is 
normalized at the AOI of the SD, and the RVS at the AOI of 
the SV is the gain difference between the AOI of the SV and 
that of the SD. Thus, the increase of the RVS at the SV does not 
necessarily mean that the gain at that AOI increases. 

Fig. 19 shows the Terra band 8 (0.412 fim) detector-averaged 
RVS derived from the EV and lunar response trending for 
MS 1 at four representative AOIs using Approach II, which 
are denoted by solid symbols. The RVS is normalized at the 
AOI of the SD. At the AOI of the SV, the RVS decreases 



Fig. 18. RVS for Aqua MODIS bands 10, 13, and 17 at the AOI of the SV. 



Fig. 19. RVS for Terra MODIS Band 8 MS 1 . Unfilled symbols: Approach I. 
Solid symbols: Approach II. 

about 10% in the first five years on orbit, then increases about 
47% in the following eight years, and exhibits a steady change 
thereafter. The change in the RVS at the SV is representative 
of the difference in the rate of gain changes derived from the 
SD and lunar observations. This observation is also evident in 
the relative RVS change presented in Fig. 16. For comparison, 
the RVS generated using Approach I is also plotted and denoted 
by unfilled symbols in the figure. The differences between the 
two sets of RVS essentially represent the differences between 
C5 and C6 since the former is based on Approach I. In general, 
a temporal increase at all AOIs is observed, with the largest 
change of around 13% observed at the AOI of the SV. The large 
RVS difference at the AOI of the SV does not imply a large 
difference in the retrieved reflectance between the two versions. 

The MODIS LIB EV reflectance or radiance is proportional 
to mi/RVS, which is inversely proportional to the gain of 
the band at any given AOI. Fig. 20 shows the corresponding 
mi /RVS trending for Terra band 8 MS 1 at the four represen- 
tative AOIs using the RVS derived from both approaches. The 
difference between the two sets of mi /RVS are expectedly neg- 
ligible at the AOI of the S V since the on-orbit gain change at the 
SV AOI is primarily determined by the lunar calibration. The 
differences vary with AOI and increase with time, particularly 
after the designed lifetime of six years. Currently, the differences 


3170 


IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 52, NO. 6, JUNE 2014 



Fig. 20. Normalized mi/RVS for Terra MODIS band 8 MS 1. Unfilled 
symbols: Approach I. Solid symbols: Approach II. 



Fig. 21. Normalized RVS for Aqua MODIS Band 8 MS 1. Unfilled symbols: 
Approach I. Solid symbols: Approach II. 

at the AOI of the nadir, SD, and end of scan of large AOI are 
about 5%, 12%, and 23%, respectively. The gain difference 
at the SD AOI indicates an increased uncertainty of the SD 
calibration for band 8 (0.412 /fin). The differences of the RVS 
at other AOIs indicate that the linear approximation for the AOI 
dependence of the RVS is not adequate to accurately describe 
the on-orbit change. Consequently, EV response trending has 
to be used to derive the time-dependent gain and RVS for Terra 
band 8. Similar behavior is also observed for other Approach II 
bands summarized in Table II. 

Fig. 21 shows the detector-averaged RVS trends from both 
approaches for Aqua band 8 MS 1 at the four representative 
AOIs. Also, as expected, there is no difference between the two 
approaches at the AOI of the SD. At the AOI of the SV, the 
difference is less than 1%, much smaller compared to that ob- 
served for Terra band 8 (0.412 /mi). Thus, the small difference 
between the approaches at the AOI of the SV indicates that 
the bidirectional-reflectance-factor [26] degradation-corrected 
SD and EV responses agree reasonably well. However, the 
differences at the AOI of the nadir and at the end of the scan 
are evident, further confirming the inadequacy of the linear 
approximation in characterizing the on-orbit RVS change. A 
similar result is also obtained for Aqua band 9 (0.443 /mi). 



Fig. 22. Terra MODIS band 8 MS 1 RVS detector differences. 

The MODIS RSB time-dependent RVS is, in principle, de- 
tector as well as band, MS, and time dependent as discussed in 
Section II. However, for most bands, the detector dependence is 
not evident, and detector-independent RVS is calculated except 
for short-wavelength bands that include Terra bands 3 and 8-12 
and Aqua bands 8-12. Fig. 22 shows the Terra band 8 MS 
1 RVS detector differences. Variations with time can be as large 
as 2.4%. For Terra band 8 MS 2 RVS, the detector differences 
are about the same as those for MS 1. For Terra bands 3 and 
9, the largest differences are 1.3% and 2.6%, respectively. For 
all other Terra and Aqua bands which have detector-dependent 
RVS, the differences are less than 1%. 

Terra and Aqua band 8 (0.412 /im) MS 1 RVS on-orbit 
changes since launch can be as large as 35% and 20%, respec- 
tively, as demonstrated in Figs. 19 and 21. For MS 2, the largest 
RVS changes are 28% and 19% for Terra and Aqua band 8, 
respectively. Apparently, Terra band 8 has a large MS difference 
while the Aqua band 8 MS difference is very small. For Terra 
band 8, the RVS on-orbit variation is even larger and can be as 
large as 47% and 38% for MSs 1 and 2, respectively, at small 
AOI since it decreased during the first few years on orbit and 
then sharply increased in the following years. For Terra bands 
9 (0.443 /im), 3 (0.469 /mi), 10 (0.488 /.mi), 11 (0.531 /mi), 
12 (0.551 /an), and 4 (0.555 /tin) MS 1, the largest on-orbit 
RVS changes are about 16%, 13%, 18%, 11%, 7%, and 4%, 
respectively, while for MS 2, the largest RVS changes for 
these bands are about 25%, 12%, 15%, 13%, 7%, and 6%, 
respectively. The RVS on-orbit variation is even larger for Terra 
band 9 MSs 1 and 2 and band 3 mirror side 2 with variations 
reaching 30%, 35%, and 28%, respectively. MS differences are 
clearly seen for Terra bands 9, 3, 10, and 4. For Aqua MODIS, 
the RVS MS difference is much smaller for all RSBs (within 
1% or less) compared to Terra MODIS RSB, and the on-orbit 
RVS change is also relatively smaller. For Aqua bands 9, 3, 10, 
1 (0.645 /im), and 2 (0.858 /im), the largest changes are about 
13%, 5%, 3%, 3%, and 3%, respectively. For all other Terra 
and Aqua RSBs, the on-orbit RVS changes are about 2% or 
within 2%. For the SWIR bands, the on-orbit RVS change is not 
noticeable, and thus, the prelaunch RVS is applied as mentioned 
in Section II. Considering that the uncertainty specification for 
MODIS RSB EV reflectance is 2%, the time-dependent RVS is 
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Fig. 23. Aqua band 3 EV reflectance trending. Solid: Time-dependent RVS. 
Unfilled: Prelaunch RVS. 

critical and essential for most of the MODIS RSB, particularly 
for the short-wavelength bands. 

Currently, six Terra bands and two Aqua bands need to use 
the desert response trending to track on-orbit RVS change while 
the other MODIS RSBs only need SD and lunar calibration 
to accurately derive their time-dependent RVS. As both in- 
struments age with passing time, either the SD calibration has 
become inadequate to catch the gain change at the AOI of the 
SD or the linear approximation for the AOI dependence for 
the on-orbit RVS change has become invalid for some of the 
bands. Consequently, Approach II has to be applied in these 
bands to derive the time-dependent RVS as previously done for 
Terra bands 1-4, 8, and 9 and Aqua bands 8 and 9. For most 
ocean bands, saturation occurs when viewing desert sites, and 
hence, ocean sites may be used to track the response change. 
However, the response trending over ocean sites are much 
noisier, and thus, it is challenging to apply Approach II to other 
ocean bands. For Terra bands 3 and 8-9, quite accurate time- 
dependent RVS has been derived from the response trending 
over desert sites as well as from the measurements of the 
onboard calibrators, but it remains a challenge to reduce the 
polarization effect at the end of the scan and the large AOI, 
particularly for forward LUT updates. 

B. EV Reflectance Trending Over Desert Sites 

The EV reflectance trending over the pseudo-invariant sites 
can be used to access the performance of the time-dependent 
RVS. Fig. 23 shows the BRDF-corrected TOA reflectance 
observed over two desert sites, Libya 2 and Libya 4, at three 
selected AOIs for MS 1 of Aqua MODIS band 3 (0.469 /xm). 
The three selected AOIs are 12.2°, 36.9°, and 48.7°, repre- 
senting the AOI of the SV, nadir, and SD, respectively. The 
reflectance trending in Fig. 23 has been normalized to the 
first measurement. The TOA reflectance trends generated using 
C6 on-orbit calibration coefficients m,\ and prelaunch RVS 
are denoted using open circles, and those generated using C6 
on-orbit calibration coefficients and on-orbit RVS are denoted 
using solid circles. Based on the assumption that the desert 
targets are temporally stable, no long-term drift is expected in 



Epoch 2000 

Fig. 24. TOA reflectance trending for Terra MODIS band 8 MS 1 . Unfilled: 
Approach I m\ and RVS applied. Solid: Approach II m\ and RVS. 

the reflectance trending, provided that the band is accurately 
calibrated. As expected, no evidence of drifts is seen at the 
AOI of 48.7° in the reflectance trending from either data sets 
since the on-orbit SD calibration is dominant in characterizing 
the gain change at this AOI. For the other two frames, the 
reflectance trending generated using prelaunch RVS exhibits 
a significant long-term drift, particularly after four years on 
orbit. The drifts accumulate to 2% to 6% at AOIs of 36.9° and 
12.2°, respectively. The above result highlights the necessity 
of a time-dependent RVS to accurately characterize the on- 
orbit sensor change at all AOIs for Aqua band 3 (0.469 /zm). 
It is demonstrated that SD/SDSM and lunar calibration can 
accurately characterize the on-orbit gain changes at the AOI 
of the SD and SV and that the linear approximation for the 
AOI dependence of the on-orbit RVS change works quite well 
for the band. A similar conclusion can also be reached for 
other Approach I bands from the reflectance trending at the 
pseudoinvariant sites using unsaturated data. 

The pseudo-invariant desert sites can be used to demonstrate 
if Approach II is necessary for deriving the time-dependent 
RVS for a MODIS RSB while the site data are not saturated. 
Fig. 24 shows the BRDF-corrected TOA reflectance observed 
over two desert sites, Libya 2 and Libya 4, at three selected 
AOIs, 12.0°, 36.6°, and 46.7°, for Terra band 8 MS 1. The cho- 
sen frames represent the SV, nadir, and SD AOI, respectively. 
The unfilled symbols denote the reflectance generated from 
Approach I, and the filled symbols denote the Approach II 
results. At the AOI of 12.0°, the observed long-term drift is 
generally within the calibration uncertainty specification of 2%. 
This also indicates that the on-orbit lunar calibration provides 
an accurate characterization of the instrument gain. The ap- 
proximately 10% long-term drift observed from the reflectance 
trending at the SD AOI further confirms the inadequacy of the 
SD to characterize the on-orbit gain change accurately. For the 
AOI of 36.6°, the unfilled symbols denote the reflectance cal- 
culated using Approach I which shows a decrease (up to 10%) 
for the first seven years followed by an increase in the trend 
thereafter. In comparison to the long-term drift observed at 
AOIs of 36.6° and 46.7°, the nonlinear dependence of the RVS 
as a function of AOI can be clearly seen. The solid symbols 
denote the reflectance trending computed using Approach II 
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Fig. 25. TOA reflectance trending for Terra MODIS band 8 MS 1. Unfilled: 
Approach I m\ and RVS applied. Solid: Approach II m\ and RVS. 

which shows a clear mitigation of the long-term drifts observed 
in Approach I. Thus, Approach II essentially is required in 
order to derive an accurate time-dependent RVS. A similar 
conclusion can be reached for Terra bands 1-4 and 9. 

Fig. 25 shows the BRDF-corrected TOA reflectance ob- 
served at the desert sites, Libya 2 and Libya 4, for MS 1 of 
Aqua band 8. Similar to Fig. 24, the unfilled symbols denote the 
reflectance generated using Approach I, and the filled symbols 
denote the reflectance produced using Approach II. At AOIs 
of 12.2° and 48.7°, the observed long-term drift is generally 
seen to be within the calibration uncertainty specification, 
which further reaffirms the ability of the onboard SD and lunar 
measurements to accurately catch the on-orbit gain change 
at the AOI of the SD and that of the SV. A long-term drift 
of about 5% is observed in the reflectance trending denoted 
by open circles at an AOI of 36.9°. This confirms that the 
linear approximation for the AOI dependence of the RVS on- 
orbit change is inadequate for Aqua band 8 and justifies the 
usage of response trending from EV targets to supplement 
the onboard measurements. The corresponding Approach II- 
derived reflectance trends show that the long-term drifts are 
significantly mitigated, and thus, Approach II is necessary for 
deriving an accurate time-dependent RVS for Aqua band 8. A 
similar result is also observed for Aqua band 9. An independent 
approach for correcting the long-term drift in Aqua MODIS 
RSB was developed by Meister et al. and a similar improvement 
was reported in bands 8 and 9 [28]. 

The EV reflectance can also be used to assess the necessity 
of applying the RVS detector dependence. Fig. 26 shows the 
differences of the EV reflectance observed by different detec- 
tors of MS 1 of Terra band 8 (0.412 /im) over the ocean site 
in MODIS C5 and C6, where the detector-independent and the 
detector-dependent RVS are applied, respectively. The TOA 
reflectance difference of each detector relative to the mean at 
the AOI of the nadir is shown in the figure. The top plot in 
the figure shows that there are clear differences among the 
reflectance observed by different detectors in Terra band 8 C5 
EV products. The differences vary with time and can be as large 
as ±2%. The bottom plot in the figure shows that the detector 
differences in Terra band 8 C6 EV reflectance products are re- 


Fig. 26. Terra band 8 C5 and C6 detector differences at the AOI of the S V. 

duced to within ±1%. A similar conclusion can also be reached 
for other bands whose RVSs are also detector dependent. 

The total uncertainty of the RVS derived using either ap- 
proaches depends on the uncertainties of input measurements 
and the associated assumptions. In addition to the detector 
dependence at short wavelengths, the uncertainty is primarily 
band, MS, and AOI dependent. Preliminary analysis shows that 
the uncertainty of the RVS derived using the first approach 
(onboard calibrators) at the AOI of the SV is generally within 
2% [24]. The uncertainty of the RVS at the AOI of the SD is 
estimated to be smaller than that at the AOI of the SV since the 
RVS is normalized at the AOI of the SD. The uncertainty of 
the RVS derived using the second approach is estimated to be 
greater than 2%, particularly at large AOI due to the impact of 
the polarization [24]. In the case of Terra band 8, the uncertainty 
of the RVS varies between 2% and 4%. A comprehensive 
uncertainty analysis of the MODIS RSB calibration will be 
discussed in detail elsewhere. 

V. Conclusion 

Some RSBs of both Terra and Aqua MODIS instruments 
have experienced significant degradation after more than 
10 years on orbit. The magnitude of the on-orbit degradation is 
more predominant in the case of Terra MODIS RSB, primarily 
due to the smoking incident prior to launch and the increased 
SD degradation since July 2, 2003. Prelaunch characterization 
exhibited a wavelength dependence in the RVS; with the largest 
change observed at the shorter wavelengths. The observed 
degradation differences at different AOIs over time are cor- 
rected using time-dependent RVS LUTs in the LIB calibration. 
Two approaches have been developed and applied to calculate 
the C6 time-dependent on-orbit RVS change. The first approach 
tracks the on-orbit RVS change using the SD and lunar ob- 
servations and EV MS ratios of the instrument response. This 
approach has been applied to derive the time-dependent RVS 
for Terra bands 10-19 and Aqua bands 1-4 and 10-19. For most 
of these bands, the RVS on-orbit change is within 5% with the 
exception of Terra bands 10-12 where it is seen to be as large as 
20%. Monitoring of the pseudoinvariant desert targets indicates 
a reasonable multiyear performance with no observable long- 
term drift. A second approach is developed to update the 
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on-orbit RVS change by supplementing the lunar measurements 
with the EV response trending from pseudoinvariant desert 
targets. This approach is used to track the on-orbit RVS change 
for the selected short-wavelength RSBs, Terra bands 1-4, 8, 
and 9 and Aqua bands 8 and 9. The RVS for these bands have 
exhibited a significant on-orbit change, about 35% for Terra 
band 8 MS 1 and 20% for Aqua band 8 MS 1 . In the case of 
Terra band 8, the difference between the two approaches of the 
RVS derivation exhibits a difference of 13% at the AOI of the 
SV. In addition to the wavelength dependence, the on-orbit RVS 
change for Terra MODIS has also indicated a MS and a detector 
dependence primarily in the short-wavelength bands. For Aqua 
MODIS RVS, both the MS and detector differences are seen 
to be comparatively within 2% for band 8 and generally less 
than 1% for other RSBs. An accurate characterization of the 
time-dependent RVS is vital in maintaining the accuracy of the 
MODIS LIB products. The lessons learned from the MODIS 
on-orbit RVS behavior can serve as a useful reference for future 
Earth-observing sensors. 
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